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Abstract 



The chiral structure of supersymmetric particle couplings involving third gener- 
ation Standard Model fermions depends on left-right squark and slepton mixings as 
well as gaugino-higgsino mixings. The shapes and intercorrelations of invariant mass 
distributions of a first or second generation lepton with bottoms and taus arising 
from adjacent branches of SUSY cascade decays are shown to be a sensitive probe of 
this chiral structure. All possible cascade decays that can give rise to such correla- 
tions within the MSSM are considered. For bottom-lepton correlations the distinctive 
structure of the invariant mass distributions distinguishes between decays originating 
from stop or sbottom squarks through either an intermediate chargino or neutralino. 
For decay through a chargino the spins of the stop and chargino are established by 
the form of the distribution. When the bottom charge is signed through soft muon 
tagging, the structure of the same-sign and opposite-sign invariant mass distribu- 
tions depends on a set function of left-right and gaugino-higgsino mixings, as well 
as establishes the spins of all the superpartners in the sequential two-body cascade 
decay. Tau-lepton and tau-tau invariant mass distributions arising from MSSM cas- 
cade decays are likewise systematically considered with particular attention to their 
dependence on tau polarization. All possible tau-lepton and tau-tau distributions 
are plotted using a semi-analytic model for hadronic one-prong taus. Algorithms for 
fitting tau-tau and tau-lepton distributions to data are suggested. 



1 Introduction 



As the LHC prepares to uncover the physics responsible for electroweak symmetry breaking, 
third-generation fermions present especially interesting possibilities as their large Yukawa 
couplings make them uniquely sensitive to details of chiral physics at the electroweak scale. 
Moreover, in the context of SUSY, there are several reasons to anticipate that SUSY signals 
may be substantially third-generation enriched. The LEP2 bounds on the masses of MSSM 
Higgs bosons suggest that cos 2/3 ~ 1, to maximize the tree-level contribution to the 
Higgs mass, favoring moderate or larger tan/3. At large tan/3, the enhanced Yukawa 
couplings of the b and r contribute negatively to the running of the sbottom and stau soft 
masses, so that sbottoms, staus, and stops all tend to be lighter than the other sfermions. 
Large Yukawas also enhance the coupling of sfermions to higgsinos, thereby increasing the 
sensitivity of sparticle decays to gaugino-higgsino mixing. In addition, large tan /3 can lead 
to large left-right sbottom and stau mixing, further lowering sbottom and stau masses [2]. 
Light stops and large sfermion mixing are also motivated by the desire to minimize the 
tuning in the quantum corrections to the Higgs mass [3J. Light third-generation sfermions 
lead to a significant enrichment of third-generation final states in the signal [2j. 

Third-generation fermions have complicated in-detector decays which require more ef- 
fort to identify and understand. On the other hand, these complicated decays allow for 
measurement of interesting properties such as polarization (tops, taus) and charge (bot- 
toms) which are inaccessible for the lighter fermions. Exploiting these properties allows for 
detailed measurement of various aspects of the Lorentz and chiral structure of new physics. 

In this paper we consider SUSY cascade decays which produce third-generation fermions 
in the final state. Our interest here is to survey the space of possibilities. Decades of work on 
supersymmetry breaking and its mediation to the Standard Model have led to a dizzying 
variety of predictions for possible sparticle spectra. At the same time, the continuing 
absence of deviations from the Standard Model predictions gives no hints as to which, if 
any, of these spectra may be preferred. In this paper we enumerate all possible patterns 
of invariant mass distributions for b-i, 1-t, and r-r pairs which arise from adjacent legs of 
on-shell SUSY cascade decays, and demonstrate how the special properties of 6's and rs 
can be used to obtain further information about the chiral properties of the MSSM. Our 
analysis is model-independent insofar as it is independent of the spectrum of superpartner 
masses; we assume only the existence of on-shell decay modes, the Lorentz structure of 
the SUSY vertices, and the field content of the MSSM [J. This type of model-independent 

^We make in addition a few mild theoretical assumptions, namely: we neglect the electron and muon 
Yukawa couplings; we neglect SUSY flavor-violating processes; we assume lepton flavor universality holds 
for electrons and muons. It is interesting to explore what happens to the intercorrelations among invariant 
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analysis has been applied to related cascade decay signals in [5]. 

Two-step on-shell SUSY cascades can give three different fundamental shapes for in- 
variant mass distributions [6J. Intermediate scalars, as in the decay chains 

X° - i^i^ - i^i^x] (1) 

give a triangular distribution for the invariant mass of the two visible standard model 
fermions, 

Here we have defined the rescaled variable 



X 



^ff,max 

Intermediate fermions, as in the decay chains 

bL,R^bx^^ ^bflL,R,biiR,L (3) 

yield "humps" 

-'^ = 4x{l-x^) = H{x) (4) 

J- (Jjtlj 

and "half-cusps" 

if ^4.'. CM (5) 

for processes without and with a helicity flip on the intermediate fermion propagator, 
respectively. These are the only possibilities when the intermediate particle is on-shell and 
the two Standard Model fermions to be combined together are adjacent in the decay chain 

We categorize the b-£ distributions which can arise from adjacent legs of SUSY cascade 
decays. In the MSSM, b-i final states can be produced from either the decay of a stop 
through a chargino or a sbottom through a neutralino to sleptons. The patterns of b-i 
invariant mass distributions arising from these cascade decays have a distinctive structure 
which discriminates between stop and sbottom initial states and establishes the spin and 
Dirac nature of the chargino. Signing muonically-decaying 6-quarks using the associated 
soft muon reveals an additional layer of structure which serves to establish the spin and 
Majorana nature of the neutralino, as well as the relative handedness of the sbottom and 
the slepton participating in the decay chain. The ability to establish spins using signed 
(for neutralino decay chains) or unsigned (for chargino decay chains) b-i distributions are 



mass distributions from cascade decays when these mild assumptions are relaxed [4], but that is beyond 
the scope of the present work. 
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entirely independent of the existence of a production asymmetry favoring squarks over 
anti-squarks. The shapes and intercorrelations of the b-i invariant mass distributions are 
sensitive to both left-right squark mixing and gaugino-higgsino mixing. 

We perform a similar categorization of the 1-t and r-r invariant mass distributions 
which can be realized in SUSY two-step cascade decays. As with b-l distributions, the 
shapes and intercorrelations of 1-t and r-r distributions are sensitive to both left-right 
squark mixing and gaugino-higgsino mixing. Invariant mass distributions involving rs 
are complicated, however, by the missing four-momentum of the neutrino coming from 
the decay of the r. The observable invariant mass distributions constructed from the r's 
visible hadronic daughters differ significantly from the underlying triangle, hump, and half- 
cusp distributions. Ditau distributions can nevertheless be used to measure superpartner 
masses with reasonable precision [71 [8]. Using the shape of r-r and 1-t invariant mass 
distributions to further establish more detailed properties of the superpartners requires the 
r visible daughter energy spectrum to be carefully taken into account. 

As the distributions of visible daughter energy depend sensitively on the polarization 
of the parent r [HI [IDl [HI [12], the observable 1-t and r-r invariant mass distributions are 
likewise dependent on the r polarization. Exploiting the dependence of the visible hadronic 
daughters on the parent r polarization enables a direct probe of the chiral structure of the r 
production vertex, opening interesting possibilities both in Higgs physics [HI [121 [131 [IH [15] 
and in the MSSM [TSl [IB [IHl [121 [20] • Careful study of 1-t and r-r invariant mass distribu- 
tions arising from cascade decays thereby opens the exciting prospect of directly probing 
the mixings and electroweak quantum numbers of the superpartners participating in the 
decayflBl [IS] • We use a semi-analytic approximation to the visible daughter energy spectra 
for the hadronic one-prong decay mode, which provides the most sensitive polarimeter for 
invariant mass distributions. Using these spectra, we plot the possible theoretical 1-t and 
r-r distributions that can arise from SUSY cascade decays, and propose algorithms for 
fitting these distributions to experimental data. 

The outline of the paper is as follows. We begin in Section [21 by considering b-l final 
states. Section [3] discusses the possible ditau and lepton-tau distributions in the limit of no 
mixing. In Section H] we turn on mixing, and discuss fitting algorithms in subsection 14.11 
Our conclusions can be found in Section [5l and details of our treatment of r decay can be 
found in the Appendix. 
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2 b-i distributions 



6-quarks are useful tools to study SUSY cascade decays. The presence of 6-tags in an event 
can help serve to separate signal from background, while the ability to sign semimuonic 
6's using the associated muon goes further to help to shed light on the Lorentz properties 
of the supersymmetric particles, as we will demonstrate. We will study the invariant mass 
distributions of b-i pairs which arise from adjacent steps in a cascade decay. 

SUSY decay chains which yield adjacent 6-quarks and leptons are the decay of sbottoms 
through a neutralino or stops through a chargino. Sbottom decays through a neutralino. 



yield both "opposite-sign" and "same-sign" (b'^^^^-i'^) final states. In the limit 

of zero squark and neutralino mixing, these processes (El) contribute either opposite-sign 
humps and same-sign cusps if the parent sbottom and the final slepton have the same 
handedness, or opposite-sign cusps and same-sign humps if the parent sbottom and the 
final slepton have the opposite handedness.. Stop decays through a chargino. 



contribute to opposite-sign final states only with no contribution in the same-sign channel. 
In the limit of zero squark and neutralino mixing, the distribution in the opposite-sign 
channel from the process (171) is a hump. 

In general events with the decay chain(s) (^U^ will also contain additional leptons com- 
ing from the subsequent decay of the sleptons, so there is some ambiguity in selecting which 
lepton to pair with the b jet. Some possible approaches to minimizing this combinatoric 
confusion have been discussed in previous studies of the related decay chain q qXi — > qil 
and its UED counterpart |6l [21], |22l, [5]. At present we will concentrate on characterizing 
the theoretical distribution for the correct b-l pairing, leaving the question of combinatorics 
to future work. 

The possible SUSY b-l distributions in the limit of no squark or neutralino mixing are 
summarized in Table [TiR If (e.g.) the decay chain bf; — >■ fe^Xi ~^ b'^Plfi exists in the signal, 
then (for a standard Majorana neutralino) so must the decay chain 6^ b^x^ ~^ b^l^lf^. 
Therefore there must exist both a hump distribution in the opposite-sign channel and a 
cusp distribution in the same-sign channel, with equal normalizations and endpoints. Since 
the hump and the half-cusp sum to a triangle, if the opposite-sign and same-sign channels 

^In the name of generality, we remark that interchanging the role of the squark and slepton does not alter 
the shapes and correlations of the h-l distributions, so that (e.g.) the decay chain b^x'i ~^ ^^'^^l 

yields the same b-l invariant mass distribution as the flipped decay chain l^Xi ^ I'^h^lr^. 



b ^ 6x- ^ ha, 



(6) 



t bxf biu, 



(7) 
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Hump Half-Cusp 

Opposite-Sign 6^ — >• 6^%° 6^ b^Xi 
it ^ b^xt 

Same-Sign 6± ^ 6±x? ^ ?^^X° 



Table 1: Possible invariant mass distributions for the h-l final states. Here by a slight abuse 
of notation ± distinguishes between (s)quark and anti-(s)quark, so that 6^ has charge ±1/3. 
Distributions for final state 1^ are obtained by exchanging hump and half-cusp. Squark 
right-left mixing is neglected and neutralinos and charginos are taken to be pure gaugino. 

cannot be distinguished, then the information about the spin of the neutralino is lost. Thus 
the ability to sign the 6-jet allows for determination both of the spin of the neutralino and 
of its Majorana nature. 

When the sbottom and the slepton have the same handedness then opposite sign dis- 
tributions are humps and same-sign distributions are half-cusps; the situation is reversed 
when the sbottom and the slepton have different handedness. It is possible to obtain 
spectra which allow all four possible decay chains bllL^R and hL^R — > hllR^i) to 

be realized simultaneously. This will complicate spin measurements as the sums of the 
overlapping distributions coming from these processes will tend to wash out the spin cor- 
relations. The greater the mass splitting between right- and left-handed squarks and/or 
sleptons, the more distinct the endpoints of the different distributions will be, and the 
easier it will be to disentangle the contributions from different processes to the total h-l 
distributions. The possibility of making a spin measurement in this channel will also de- 
pend on the relative branching fractions into sleptons of different handedness. Consider 
the case when right- and left-handed sleptons are nearly degenerate. Then the differ- 
ent gauge quantum numbers of the right- and left-handed sleptons will still yield (e.g.) 
T(hL bXi hllL)/T(hL bXi blln) ^ 1. The total h-l distribution will then retain 
some spin information. 
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Stops, on the other hand, contribute only to the opposite-sign channel, so comparison of 
opposite-sign and same-sign distributions will be an important tool to disentangle possible 
contributions from stop and sbottom squarks. 

When the intermediate chargino is pure gaugino, only decay chains beginning from 
initial ti can contribute to b-l final states. However, as the lepton Yukawas are negligible, 
it is also sensible to consider the case where the chargino is nearly pure higgsino but 
still decays to lepton-slepton pairs as a gaugino. In this case, the can decay through 
an up-type higgsino in the decay chain b'^xf ~^ If^l^V^ yielding an opposite-sign 

hump distribution. The can decay through a down- type higgsino via the decay chain 
if^ — *• b'^xt ~^ b^l^V, which yields an opposite-sign half-cusp distribution. 

In stop decay, the 6-tag alone suffices to establish that the intermediate fermion is 
Dirac, rather than Majorana, and the additional sign information is not necessary: if the 
intermediate fermion were Majorana, then summing the distributions for the final states 
b-l^ and b-l^ would give a triangle, while for an intermediate Dirac particle, only one of the 
hehcity states can contribute, leading to a hump distribution (or a half-cusp distribution 
when the chargino is a down-type higgsino) even when both final states are summed. 

This is an example of a general point: it is possible to observe nontrivial (non-triangular) 
distributions in the absence of sign information, if a symmetry forbids one channel from 
contributing (as for stop decays through charginos, here, or as for Dirac neutralinos [1]). 
Nontrivial distributions also can be obtained if the two channels contribute with unequal 
weights, as is the case when a production asymmetry favors squarks over anti-squarks [6]. 
However, event-by-event signing of the semimuonic 6's allows a direct observation of any 
angular correlations coming from intermediate Majorana neutralinos, independent of any 
possible production asymmetry. 6-jet signing as a tool to improve spin measurements at 
the LHC has been mentioned in [21] . 

We now go on to discuss how this story is modified in the presence of nontrivial squark 
and neutralino mixing. 

2.1 h-l distributions with nontrivial squark and neutralino mixing 

Left-right sfermion mixing renders the SUSY sfermion-fermion-gaugino vertices less chiral 
and thereby alters the observable kinematical distributions [25]. Similarly, higgsino inter- 
actions proportional to a fermion Yukawa coupling involve the opposite chirality of the 
fermion relative to gaugino interactions. Sfermion and neutralino mixings, then, have a 
qualitatively similar effect on the invariant mass distributions. In this section we will de- 
tail the sensitivity of difermion invariant mass distributions to both sfermion mixing and 
neutralino mixing, including both effects simultaneously. While here we concentrate on 
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6-quarks and b-l distributions, the same physics will be relevant to lepton-tau and ditau 
distributions in section HI 

Consider the hump and half-cusp distributions which arise from decay chains with an 
intermediate fermion. Once mixing is turned on, both helicity states of the intermediate 
fermion can contribute to a given channel, with relative weights determined by the mixing. 
The observable same-sign and opposite-sign distributions are then a weighted sum of hump 
and half-cusp distributions. Define the squark left-right mixing angles through 

= cos e~q*ji + sin e~qL. (8) 

(As usual, qi is taken to be the lighter of the two squarks. In our conventions is a 
left-handed anti-quark, and thus is an anti-squark.) Define also the unitary matrix U 
which diagonalizes the neutralino mass matrix. The sbottom-bottom-neutralino couplings 
are then governed by the interaction Lagrangian 

C = h + ix'^M J + h {bnx'^ + {x'^H n^J + H.c. (9) 

with parameters 

= sin 9iX,U:, + cos e^y^U*^, (10) 

< = sine-,(^-^U,w + ^U,B^+cose-,X,U*, (11) 

n£ = cose-,X,U*,~sinej,^U*R„ (12) 
n^i = cos 01 ^Uiw + ^^UiB^ -sin e^XbU^i (13) 

Note that n^- — > 0, raf ^ as both mixings are turned off, that is, as 9^, Udi 0, while 
j, nfj remain finite. The index i specifies the neutralino mass eigenstate i, while the 
indices d, B, W run over the gauge eigenstates (here down-type higgsino, bino, and wino, 
respectively). We now define the angles 
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COS ai, = , o I, ,r I, , cos a2i = , . „ , ,0 ; (14) 



in the limit of zero mixing, cos an and cos a2i both go to unity. Decays of sbottoms through 
the neutralino will be weighted by these angles. Meanwhile, the stop-bottom-chargino 
couplings are governed by the interaction Lagrangian 

C = h {bRX- 4 + ixtVbl cl) + h {hRX- 4 + ^xt)K 4.) + H.c. (15) 
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with parameters 



i = - sin 9tgViw + cos OiXtViu (16) 

= smOfX^W:, (17) 

C2,i = - COS OigViw -sin OiXtViu (18) 

c£ = cos 9iX,W^„ (19) 



The unitary matrix V (not the CKM matrix!) diagonahzes the positively-charged left- 
handed charginos {W^, h^), and the unitary matrix W diagonahzes the negatively-charged 
left-handed charginos {W~, h^). We then define the angles 

I ■ I ^ I ■ I ^ 

COS^A, = I B|2 4.'. L 12 ' COS^/^2^ = . r '.^ H 12 ' (^0) 

where again cos and cos /?2j go to unity in the limit of zero mixing. Decays of stops 
through the chargino xt be weighted by these angles. Table [2] summarizes the 6-lepton 
distributions in the presence of mixing; a similar approach to mixing was taken in |5]. 

Notice that for stop squark decays through charginos, admixture of the cusp distri- 
bution depends on a sizable down-type higgsino component of the intermediate chargino. 
This can be seen as follows. As the lepton Yukawas are negligible, the chargino decay to 
lepton-sneutrino pairs will occur through its gaugino component. To obtain a half-cusp 
distribution, the b must then be right-handed (as its charge is fixed). But as the wino 
does not couple to the right-handed b, a half-cusp distribution can only be obtained if the 
intermediate chargino has a down-type higgsino component, independent of any possible 
mixing of the stop squarks. The coupling of the higgsino to the bji depends on and is 
therefore enhanced at large tan/3. In short, if the chargino is taken to have no down-type 
higgsino component, then stop mixing alone will not alter the observed hump distribution. 
In this context it is useful to recall that no sign information of the 6-quark is necessary to 
observe the hump distribution coming from stop decay. 

Fitting the observed b-i distributions to sums of humps and half-cusps thereby measures 
the mixing parameters of equations f|T^ and fl20|) . These mixing parameters are the only 
linear combinations of Lagrangian parameters which can be measured using these cascade 
decays alone. For stop squark decays, the mixing parameters and any production asym- 
metry can be separately measured by considering the separate distributions b-l^ and b-l^ . 
For sbottom squark decays, the mixing parameters and the production asymmetry can 
be separately measured if the 6-quark is signed and the distributions for all four possible 
combinations of lepton and b signs are independently examined. 
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Hump Half- Cusp 

Process Weight Process Weight 



^ b^l^l^ sin^ an ^ b'^l^l^ cos^ an 

^ b^l^l%^ cos^ an ^ b^l^l^ sin^ an 

Opposite-Sign 6f b'^Xi bf b'^x'i 

^ b'^Plf cos^ a2i ^ b^l^lf^ sin^ a2i 

b^l^l^ sin^ a2i b^l^l^ cos^ 



6=F/±z/L cos^/5(i,2)j ^b'^l^VL sin^/5(i,2)i 



•^b^l^ll cos^aii ^b^l^ll sin^aii 

b^^l l]i sin^ > b^l^V^ cos^ an 

Same-Sign 6^ ^ ^ 

b^l l[ sin^ a2j b'^l^V[ cos^ 0:21 

b^l cos^ a2i ^ b^l^r^ sin^ a2i 



Table 2: Distributions for the b-l invariant mass including both gaugino-higgsino mixing 
and third-generation squark mixing. Here by a slight abuse of notation ± distinguishes 
between (s)quark and anti-(s) quark, so that b^ has charge ±1/3. The relative weights are 
normalized such that the coefficients for the hump and cusp distributions in a given channel 
sum to unity. The angles aki, Pki are defined in equations (fn|) and fl2U|) . The endpoints 
of the distributions will depend on the masses of the superpartners participating in the 
cascade. 
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2.2 6-jet signing and measurement of mixings 

In practice signing the 6-quark cannot be done with arbitrary purity, and there will be 
a non-trivial but well-characterized mis-sign rate. In particular there is an irreducible 
contribution to the mis-sign rate coming from oscillations of the parent 6-quark inside 
neutral 5-mesons, which is of order 12% p6]. Recent simulations indicate that a mis-sign 
rate of 15.4% can be achieved, with an efficiency of 1.2% [27]. To measure spin and mixings, 
b-i distributions should be fit to sums of humps and half-cusps. Let the mis-sign fraction 
be F, and consider decay of 6-squarks through neutrahnos. If the theoretical distribution 
in a given channel is 

Dtheoryix) = fH{x) + {1 - f)Cix), 

where / and (1 — /) are the (sines and cosines of) angles listed in Table [H then the 
experimentally observed distribution is 

DoUx) = {l-F){fH{x) + {l-f)C{x)) + F{{l-f)H{x) + fC{x)) 

= Mx) + (l-/)C(x), (21) 

where 

f = f + F- 2fF (22) 

Fits to data then directly measure /, which through fl22l) measures the mixing parameters 
of equation (IT^ . In a realistic situation there will be a trade-off between purity, that is, 
minimizing F, and acceptance. In a full analysis any p^-dependence of F can be included. 

In Figure [1] we plot the observable b-i distributions for an intermediate pure bino and 
cos6'f, = 0.775, assuming no production asymmetry. Two sets of curves are shown. The 
outer (red) pair are the opposite- and same-sign distributions which would be observed with 
a 15% mis-sign rate. The inner (blue) pair are the opposite- and same-sign distributions 
which would be observed with a 30% mis-sign rate. Even with a 30% mis-sign rate the 
deviation of the distributions from each other and from the triangular distribution is clear, 
indicating the presence of an intermediate Majorana fermion. Decreasing the mis-sign rate 
further to 15% significantly enhances the difference between the two channels, and hence 
the sensitivity to mixing. 

Let us finally mention that in the semi-muonic decay modes of the b, the missing energy 
carried away by the neutrino alters the energy distribution of the observed 6-jets, and a full 
analysis must account for this effect. As the polarization of the 6-quark is almost entirely 
randomized in hadronization [2B], the visible spectrum of 6-quark decay products is not 
dependent on the details of the vertex where the b originated. The net effect of the loss 
of the neutrino four-momentum can be described by a convolution of the invariant mass 
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Figure 1: Observable b-i invariant mass distributions arising from b decay through a neu- 
trahno, for cos6'f, = 0.775 and intermediate pure bino. The outer (red) pair of curves show 
the observable distributions assuming a 15% mis-sign rate. The inner (blue) pair of curves 
show the same distributions assuming a 30% mis-sign rate. 

distributions involving the parent 6-quark with a calculable universal transfer function 
which accounts for the distribution of visible energy in the &-jets after the loss of the 
neutrino. 



3 Ditau and lepton-tau distributions: no mixings 

As taus decay within the detector via r uX, the full four-momentum of the r is not 
observable, and events involving final state rs require careful attention to characterize and 
understand. However, the hadronic decay modes of the r do allow right- and left-handed rs 
to be distinguished statistically. A left-handed r preferentially emits the neutrino parallel 
to its direction of motion, resulting in a softer spectrum of visible decay products, while 
a right-handed r preferentially emits the neutrino anti-parallel to its direction of motion, 
resulting in a harder spectrum of visible decay products. This difference in the energy 
distributions of the r daughter products therefore can be used as a handle on the possible 
chiral couplings of new physics to r leptons [TOl [HI [12] . 

When the r is highly boosted in the lab frame, as it is in most cascade decays, its visible 
decay products d are coUinear, and to good approximation we can take 

Pd = ZPr, 

where z is the fraction of the lab frame r energy carried by the daughters d, 
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Figure 2: The r daughter energy fraction distributions -P(i)(-2) for hadronic one-prong taus, 
for positive (red) and negative (blue) hehcity taus. 

The quantity z is invariant under boosts along the direction of the tau. Also, when the tau 
is highly boosted, helicity and chirality may be used interchangeably. The distribution of 
z is correlated with handedness of the r, as noted above. 

At the level of the parent rs, two-step on-shell SUSY cascades can lead to triangle, 
hump, and half-cusp invariant mass distributions. Intermediate staus, as in the decay 
chains 

X° - rl^Tl^ ^ rl^Tl^x] (24) 
give triangles, while the decay chains 

lL,R^iX^,^trTi^^R.iTTn,L (25) 

and 

tl,r ^Xi tttr^l (26) 

yield humps and half-cusps. For final state r's, these distributions are not directly ob- 
servable. Rather, we must convolve these distributions with the probability P^-^iz) that a 
parent r decays to visible daughter particle(s) d with momentum fraction z = Pd/Pr- The 
probability P^-^ depends on the hehcity of the parent r, denoted ±. CP invariance ensures 
that the energy distributions of a r with a given helicity and its anti-particle are identical. 
Thus the negative-helicity and the positive-helicity r"*" have identical energy distribu- 
tions, which we denote by P^j^{z) throughout. Similarly, the energy distributions for the 
positive-helicity r~ and its antiparticle, the negative-helicity r"*", are given by P^^-^{z). 

Our analysis uses hadronic one-prong taus as they yield the greatest sensitivity to the 
tau polarization; our treatment of the one-prong decay mode closely follows that of [12j. 
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We model the energy fraction distributions P^^{z) for one-prong taus by summing the 
three dominant contributions to this mode, the decays un^, up^ wk^tt^, 

and r~ va^ utt^ti^tt^. Details of this computation are presented in the Appendix. 
The total daughter energy fraction distributions P^-^ (z) that result are shown in Figure El 

The functions P^^{z) for three-prong taus are comparatively insensitive to tau polar- 
ization. The three-prong decay mode is principally due to the decay r aiu, and the 
total visible energy fraction for the three-prong decay mode is approximately equal to the 
visible energy fraction contributed by the ai's to the one-prong decay mode. The mass 
difference between the ai and the r is not large, and the contributions of longitudinally 
and transversely polarized ai mesons add to a nearly spin-independent quantity (further 
details can be found in the Appendix). Therefore including the three-prone decay mode 
does not increase sensitivity to r polarization in invariant mass distributionqj- 

In lepton-tau final states, the observable invariant mass variable is rriM, where d again 
denotes the visible r decay product (s). The variable mj^ is distributed according to 



1 dr I dz „^ , ^ dV 



Tdmja T Jml^ z ^'^'^ ' dm} 



(27) 



=mlJz 



In ditau final states, we have similarly 

1 dV ^ dzi dz2 r,± / ^ dT 



r^m2 ~ r / 2 / 2 / y P(d)i^v ^ P(d')K^v > 2 



(28) 



Equations (127]) and fl28l) lead to calculable predictions for the invariant mass distributions of 
detected tau decay products which will in general depend on the underlying SUSY process, 
the polarization of the parent tau, and the decay mode(s) selected. We plot invariant mass 
distributions of this form normalized to unity, rather than to the 45% branching ratio (or 
the branching ratio squared) into the hadronic one-prong decay modes we model. 

Our interest in the rest of this paper will be to explore how information about r po- 
larization can be used in conjunction with invariant mass distributions to further measure 
properties of a general SUSY model. Our principal aim here is to establish the range of 
theoretical possibilities. 

For the purposes of the remainder of this section, we assume purely chiral couplings, 
that is, we work in the limit of vanishing left-right stau mixing and Yukawa couplings, and 
assume that the neutralinos participating in the cascade decays couple purely as gauginos. 



•^The three-prong decay modes can, however, be used as an effective polarimeter by examining the 
distribution of energy among the daughter pions, which is sensitive to the polarization of the ai, and 
therefore to the heUcity of the t [M] . Cuts on the relative energy distributions of the daughter hadrons as 
a way to distinguish between tau polarizations in SUSY cascades have been discussed in pO] , 
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Figure 3: Ditau triangles, for intermediate tr (red) and intermediate tl (blue), using one- 
prong hadronic rs. The distributions are normalized to unity. No mixing. 

This is an idealization, and in many well-motivated scenarios left-right stau mixing is 
appreciable. First we will consider ditau final states, then in section 13.21 move onto lepton- 
tau final states. Then in section H] we will incorporate the effects of mixing on both types 
of final states. 

3.1 Ditau distributions 

Neutralino decay through staus, 

X° - tIrtIr - rl^rl^xl (29) 

yields a triangular distribution for the ditau invariant mass in the opposite-sign channel. 
We plot the resulting observable invariant mass distributions of the tau decay products 
in figure [31 These curves are a result of convolving the underlying triangle distribution 
with (1) P~ {zi)P~ {z2)., for an intermediate fi, and (2) P^{zi)P'^{z2)., for an intermediate 
Tr. These curves are fairly well separated and it should be possible to clearly distinguish 
between these two scenarios. 

Neutralinos intermediate between two staus lead to an underlying hump distribution 
in the same-sign channel, simultaneous with an underlying half-cusp distribution in the 
opposite-sign channel. These humps and half-cusps involve one tl and one tr and therefore 
the observable distributions are obtained by convolution with the (H — ) combination of 
tau energy transfer functions. The resulting distributions are plotted in figure |H These 
distributions must have equal normalization and endpoints. While in practice the upper 
endpoint may be difficult to discern, this nonetheless translates into a stringent correlation 
on the relative locations of the peaks in the opposite-sign and same-sign distributions. 
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Figure 4: The ditau hump (green) and half-cusp (orange), convolved with the (H — ) com- 
bination of energy distribution functions. Using one-prong hadronic rs. No mixing. 

An important question is at what level the possible opposite-sign distributions can be 
distinguished from each other. Toward this end we plot both possible triangle distribu- 
tions and the half-cusp in Figure [51 While with enough statistics the three curves might 
possibly be distinguished, on a practical level generating and testing hypotheses to ex- 
plain an opposite-sign ditau signal will proceed first by cross-channel comparisons. First, 
if the opposite-sign ditau distribution is an underlying half-cusp, then there should be a 
corresponding underlying hump distribution in the same-sign channel. The absence of a 
same-sign signal strongly suggests that the opposite-sign signal is an underlying triangle 
distribution, due to intermediate staus. In addition, opposite-sign triangles and opposite- 
sign half-cusps fit into hypotheses which make different predictions for object counts in the 
rest of the signal events. In particular, if the opposite-sign signal is due to half-cusps, then 
there should be a larger number of leptons in the event (most likely r's) coming from decays 
into the r initiating the decay chain, and from subsequent decay of the f terminating the 
decay chain. These additional r's will naturally present some combinatorial complications, 
which again we will not address here. 

3.2 Lepton-tau distributions 

Lepton-tau distributions arise from slepton decay to a stau through a neutrahno, 

or the analogous process with initial stau and final slepton. In the absence of mixing, there 
are four possible observable distributions in lepton-tau channels, namely the hump con- 
volved with the (+) energy distribution function; the hump convolved with the (— ) energy 
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Figure 5: Opposite sign ditau distributions: both possible triangle distributions (with 

tau polarizations (++) in red and ( ) in blue), plotted against the half-cusp (with tau 

polarizations (H — ), in orange). Using one-prong hadronic rs. No mixing. 

distribution function; the half-cusp convolved with the (-I-) energy distribution function; 
and the half-cusp convolved with the (— ) energy distribution function. These distributions 
are plotted in Figure [61 Depending on the handedness of both the slepton and the stau 
participating in the decay chain, all of these shapes may appear in either the same-sign 
channel or the opposite-sign channel, as we summarize in Table [31 Note that, for any 
given process, there is a specific prediction for which of the four shapes must appear in the 
same-sign channel and which must appear in the opposite-sign channel. They must have 
the same normalization, and the same endpoints. In addition it is worth pointing out that 
the helicity of the tau in both the opposite-sign and same-sign channels is the same. 

It should be readily possible to distinguish whether the cusp distribution occurs in the 
same-sign or the opposite-sign channel, as even after convolution the hump and half-cusp 
distributions are fairly distinct. This allows one to distinguish between the scenario where 
the slepton and the stau have the same handedness, and the scenario where the slepton and 
the stau have opposite handedness. To proceed further one would like to identify the hand- 
edness of the stau and therefore of the slepton. This requires comparing the Cr distribution 
to the Cl distribution, and likewise between the Hr and the Hi distributions. While this 
may be challenging for the humps, the cusps are more distinct. The discriminatory power 
is enhanced by the existence of two channels which must both have the same polarization. 
With enough statistics we expect that the identity of the stau and therefore of the slepton 
can be discerned. 
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Figure 6: All four possible lepton-tau distributions: the hump with negative (green) and 
positive (orange) polarizations, the cusp with negative (blue) and positive (cyan) polariza- 
tions. Using one-prong hadronic rs. No mixing. 



Process Same-Sign Opposite-Sign 



Ir- 




Cr 


Hr 


Ir- 




Hl 


Cl 


II- 




Hr 


Cr 


h- 




Cl 


Hl 



Table 3: Possible lepton-tau distributions in the absence of mixing. Here we denote by 
Cr the half-cusp distribution convolved with the positive polarization energy distribution 
function, and by Hl the hump distribution convolved with the negative polarization energy 
distribution function. Identical conclusions pertain if the roles of the stau and the slepton 
are reversed. An experimental determination of the channel in which the half-cusp and 
hump distributions appear determines the relative handedness of the slepton and stau. 



17 



4 Ditau and lepton-tau distributions in the presence 
of mixing 

In many realistic SUSY models, mixing in the tau sector is nonnegligible. As left- and 
right-handed r's have different daughter energy spectra, the observable invariant mass 
distributions will be a weighted sum of the distributions for purely left- and purely right- 
handed rs, with weights determined by the mixing parameters. In addition, as we saw 
with b-l distributions in section [2], reducing the chirality of the fermion-sfermion-neutralino 
vertices serves to wash out angular correlations from intermediate fermions. A careful fit 
of observed tau distributions therefore has the potential to probe the chiral structure and 
mixings of the new physics. 

The Yukawa interactions between the two stau mass eigenstates, the i*^ neutralino, and 
the right and left-handed taus are 

A„t = n {trx'I y?, + y^,) + T, (r^x° + y^,) + H.C.. (30) 

(In our conventions, r^j is a left-handed anti-tau.) The stau-tau- neutralino Yukawa cou- 
plings are 

sin e~XrU^i + cos e~V2g'U*B, (31) 



yl,^ = sine- [-:^U,w - -^UiBj + cos^^^A.f/d. (32) 

= cos^~A.t/,*,-sin^~v^(7'f/B. (33) 
y^. = cosO- (^--^Uiw - ^UiB^ -sin e-XrUdi , (34) 

with the property that ?/2i ~^ l^l/fi ~^ l^^/fi ~^ O^Z/f^i as the both the stau and 
neutralino mixings are turned off, that is, as 6~, Udi — > 0. The combinations of these 
parameters which enter into the observable distributions are the relative probabilities of 
producing right- and left-handed taus at each vertex, 

COsVl,i = I R 12 Y I L 12 (35) 





\yU 


2 


\yU 


2 + 


\y2^ 





COS 02,^ = ,r 10 .P|o • (36) 



4.1 Ditau triangles 



Consider first the FSF process, neutralino to stau to neutralino. This gives a ditau "tri- 
angle" . It is possible now for each r to be either positively or negatively polarized, with a 
probability depending on the mixings in both the neutralino and the stau sectors. 
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Figure 7: All three possible ditau triangles: Tll (blue), Tlr = Trl (green), and Trr (red). 
Using one-prong hadronic rs. All curves have the same normalization and endpoints. 

The general distribution for this process is a sum of three possible fundamental shapes, 
corresponding to the triangle distribution convolved with the three different combinations 
of tau polarizations, Tll.Tlr.Trr. These shapes are plotted in Figure [71 For general 
mixings, the cascade t^i't '^'^X^ iiow gives the invariant mass distribution 



This defines a two-parameter family of distributions. In principle one could fit the observed 
ditau invariant mass distribution to this formula, using a two-parameter fit. This would 
measure the mixings in the stau and neutralino sectors as well as the helicities of the taus. 
Such an analysis is difficult, however, as the fit is not completely straightforward. First, 
the distribution T^r only differs from the average (Tll + Trr)/2 by a few percent. Second, 
it is difficult in practice to accurately locate the upper endpoint of these distributions. We 
now discuss a method of analysis which deals with these two issues. 

The ditau distributions are weighted towards smaller invariant mass, x < 0.5, even 
for underlying half-cusp or triangle distributions which peak at large invariant mass, due 
to the energy lost to neutrinos. An unfortunate consequence of this feature is that the 
upper endpoints are poorly defined for many of these distributions, as can be seen in 
Figures [3, [HI and (TDl Anchoring the distributions by their endpoints may then not be 
feasible experimentally. The most visible feature of all of these distributions is not the 
upper endpoint but rather the location of the maximum. Moreover, the shape and relative 
location of the maxima are distinct features of the various distributions. A reliable way to 




(37) 
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fit experimental data to (sums of) tliese distributions is tlius to fit to tlie location of tlie 
peak, rescaling tlie normalizations to preserve the total area of the distributions — that is, 
the total number of events — while allowing the relative locations of the endpoints to vary. 
Given a measured ditau distribution with a peak located at invariant mass 

mkh' = nT-peak, (38) 

the experimental data can be compared to theoretical distributions rescaled to have the 
same peak location. 

To do this, consider a generic ditau distribution T{a}{x) which is a linear combination 
of the basic theoretical distributions Tll, Trr and Tlr, 

T{a}{.x) = aLLTLL{x) + ajiRTnR{x) + aiRTLnix). 

This distribution T{a}(x) has a maximum at a fixed numerical value Xpeak = C{a}, which is 
uniquely determined by the coefficients aiL, (Irr and aiR, Recall that the scaled variable 
X is the ratio of the invariant mass to the endpoint, 

mhh' 
X = . 

^end 

If the (scaled) theoretical distribution T{a}(x) is to describe the data with a peak at mp^ak, 
it must have an endpoint at niend = 'f^peak / C{a} ■ Define a new variable y which is scaled by 
the visible location of the peak, rather than the location of the endpoint: 

y ^ J>1_ _ mhh' _ mhh' 

C{a} C^a}miend mtp^af^ 

In terms of y, the properly normalized theoretical distributions to fit to data are 

{c{a}y) , (39) 

where /q^^'^^"* dyP^ajiu) = 1- By construction all distributions P^a} have a peak at the same 
location. The three rescaled triangle distributions are plotted in Figure [HI 

This rescaling procedure maximizes the distinguishability of the different possible dis- 
tributions and allows the fit to be performed without any knowledge about the location of 
the endpoint in the experimentally observed distributions. However, it does require back- 
grounds to be well-characterized, as the total number of events needs to be well understood. 

We now address the approximate degeneracy between T^n and (Tll + Trr)/2. Writing 
sin^0i^j = Wj, equation fl371) can be rewritten 

D{X) = Trr{x) + [Wi + W2)(Tlr - Trr) + WiW2{Trr + Tll - 2Tlr). (40) 
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Figure 8: The three possible ditau triangles: Tll (blue), Tlr = Trl (green), and Trr (red), 
rescaled to have the same normalization and location of peak. Using one-prong hadronic 
rs. 

The final term in parentheses can be neglected, subject to the limits of experimental pre- 
cision. Moreover, this term is multiplied by a coefficient quadratic in Wj. If the angles 0ij 
are small — so that the r's being produced are predominantly right-handed — then the final 
term in equation (HUj) is doubly small. The experimental data can then be fit to 



using a one-parameter fit. This is a useful parameterisation when one suspects that the 
taus being produced are mostly positively polarized, as one would be able to learn from 
an examination of the relative energy carried by charged and neutral hadrons in the recon- 
structed taus. If the taus are mostly left-handed then a more useful parameterization of 
equation ^S7\ is obtained by taking cos^ 0ij = Wj, and fitting to the resulting approximate 
distribution 



Of course, in an intermediate situation the parameterization cos^0i,i = tfi, cos^ 01,2 = 
1 — W2 may be more convenient. In all of these cases the fit measures the sum of the 
relative probabilities of producing right-handed versus left-handed taus at both vertices. 

4.2 Ditau humps and half-cusps in the presence of mixing 

In the presence of mixing, the cascade 




(41) 




(42) 




TTTi 



(43) 
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Figure 9: The three distributions which contribute to the same-sign ditau distribution: 
Hlr (black), Cll (blue), and Crr (red). Using one-prong hadronic rs. All curves have the 
same normalization and endpoints. 

now gives the invariant mass distribution 

Dssix) = (cos^0i_icos^02,i +sin^0i,iSin2 02,j)^^LR(a;) (44) 
+ cos^ sin^ 02,j Crr + sin^ cos^ (f>2,i Cll{x) 

in the same-sign channel, and the distribution 

Dos{.x) = (cos^0i,icos^02,i + sin^0i,iSin^(/)2,j)CLfi(x) (45) 
+ cos^ sin^ 02,i Hrr + sin^ cos^ 02,j Hll{x) 

in the opposite-sign channel. As we have discussed above, object counts and cross-channel 
correlations will likely first be used to distinguish between FSF scenarios and SFS scenarios, 
so for the moment we concentrate on the mixing-induced modifications to the ditau distri- 
butions which result from the process fl43p alone. The three distributions which appear in 
equation fH4|) are plotted in Figure M The distribution Hlr, which would appear in the 
limit of no mixing, is shown in black, and the mixing-induced contributions Cll and Crr 
are shown in blue and red respectively. Similarly, the three distributions which appear in 
equation are plotted in Figure [TDl The distribution Clr is shown in black, and the 
mixing-induced contributions Hll and Hrr in blue and red respectively. 

To compare the distributions to data we rescale the distributions to fit the location of 
the peak and the total number of events, as with the triangles. The rescaled distributions 
are plotted in Figures [11] and [121 Unlike the triangles, all three curves in each channel are 
distinct enough to allow the possibility of a full two-parameter fit, although high statistics 
would be required. In addition one is able to perform the same measurement in both the 
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Figure 10: The three distributions which contribute to the opposite-sign ditau distribution: 
Clr (black), Hn (blue), and Hrr (red). Using one-prong hadronic rs. All curves have 
the same normalization and endpoints. 



opposite-sign and same-sign channel. In this case the probabilities of producing right- and 
left-handed taus at each vertex could separately be measured, in principle. In practice, any 
measurement of this sort would require a good strategy for dealing with the combinatorial 
complications of additional rs in the event, as well as the likely presence of opposite-sign 
ditau triangle distributions in addition to the opposite-sign (mixed) half-cusp of fj45|) . 

4.3 Lepton-tau distributions with mixing 

The possible lepton-tau distributions in the absence of mixing are listed in Table [31 Once 
mixing is turned on, a channel which began with a hump (half-cusp) distribution and a 
given polarization in the absence of mixing will also have a contribution of the half-cusp 
(hump) distribution, convolved with the opposite polarization. Thus, for instance, the 
process 



As there is only one r vertex, this is only a one-parameter fit. The result of the fit is the 
relative probability of emitting right- versus left-handed taus at the tau-stau-neutralino 



/rri 



1 



gives the distribution 




in the opposite-sign channel, and the distribution 
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Figure 11: The three distributions which contribute to the opposite-sign ditau distribution: 
Clr (black), Hll (blue), and Hjm (red), rescaled to have the same normalization and 
location of peak. Using one-prong hadronic rs. 
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Figure 12: The three distributions which contribute to the same-sign ditau distribution: 
(black), Cll (blue), and Crr (red), rescaled to have the same normalization and 
location of peak. Using one-prong hadronic rs. 
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Figure 13: The distributions Hr{x) (red) and Cl{x) (blue). Using one-prong hadronic rs. 

vertex; the opposite-sign and same-sign channels independently constrain this probability. 
Notice that the lepton-tau distributions appear in pairs: in any given channel, the distri- 
bution must be a weighted sum either of H^{x) and Cl{x) or of Hl{x) and Cr{x). There 
remains the discrete ambiguity of which of these pairs of distributions occurs in the same- 
sign and which in the opposite-sign channel. However, as the lepton-tau distributions are 
more distinct from one another, and as less energy is lost to neutrinos, one scenario should 
be clearly preferred. 

In Figures [12] and [H] we plot the two separate combinations of lepton-tau distributions 
which can appear in any given channel. Compared to the ditau distributions, the peaks of 
these distributions are broader, and the endpoints of these distributions are more distinct. 
In these channels, fits which anchor the endpoint of the distribution rather than the peak 
may yield better results. For comparison we plot the rescaled distributions with identical 
peak locations in Figures [15] and [THl 

5 Conclusions 

We have categorized the possible h-l, t-1, and r-r invariant mass distributions which can 
arise from two-step on-shell SUSY cascade decays and outlined several ways in which the 
special properties of 6's and r's can be used to probe detailed properties of the underlying 
MSSM Lagrangian. 

In the MSSM, b-l distributions arise from the decay of a squark through a neutralino 
or chargino to a slepton. The interplay between the nontrivial spin of the intermediate 
fermion and the chiral vertices leads to a rich structure of invariant mass distributions. The 
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Figure 14: The distributions Cr{x) (red) and Hl{x) (blue). Using one-prong hadronic rs. 




Figure 16: The distributions Cr{x) (red) and Hl{x) (blue), rescaled to have the same 
normahzation and location of peak. Using one-prong hadronic rs. 

distinctive structure of the invariant mass distribution arising from stop decay through a 
chargino establishes the spin of the chargino and depends on the mixings in the chargino- 
stop sector. Signing the 6-jet further allows the separate contributions of the two helicity 
states of a Majorana neutralino to be resolved, and opens the door to measurements of the 
neutralino spin and mixings in the sbottom-neutralino sector. 

Invariant mass distributions for final state rs depend on the polarization of the r. For 
analyses using invariant mass distributions only, hadronic one-prong rs provide the greatest 
sensitivity to r polarization. We have used a semi-analytic model of the one-prong decay 
mode incorporating r decays to 7r,p, and ai, which together account for 90% of all one- 
prong r's. With this model we computed the possible theoretical distributions for the 
observable t-1 and r-r invariant masses. 

The ditau triangle distribution has excellent prospects for yielding detailed measure- 
ments under realistic conditions, due in part to the minimal combinatorial confusion in 
such events. Precision fits of the ditau triangle can directly establish the handedness of 
the intermediate stau and measure mixings in the stau-neutralino sector. This fit is not 
completely straightforward as, due to the two missing neutrinos, the upper endpoint of the 
ditau triangle distribution is difficult to discern. The peak of the ditau triangle is more 
readily located, and we have proposed an algorithm which fits ditau distributions to the 
location of the peak, rather than the location of the endpoint. 

The distributions which we have examined here are uniquely sensitive to spins and 
mixings and present several interesting possibilities for measurements. Much more work 
is necessary, however, in order to assess when and how well such measurements can ac- 
tually be made in realistic circumstances. One of the chief obstacles to fully utilizing the 



27 



correlations identified here is the combinatorial confusion associated with identifying the 
correct pairing of objects in an event. Without a good strategy to minimize this ambiguity, 
the ability of (e.g.) b-l invariant mass distributions to determine neutralino spin is greatly 
reduced. As an acceptance price must be paid for signing the 6-jets, cuts to reduce combi- 
natorial background must be applied with care in order not to run out of signal. However, 
there are several possible approaches to reducing the combinatorial background and we are 
optimistic about the possibilities of performing interesting spin and mixing measurements 
under realistic conditions. We intend to return to this point in future work. 

In our systematic consideration of b-l, t-1, and r-r invariant mass distributions, we 
have mapped out the space of possibilities for what shapes and intercorrelations of these 
distributions can appear, and for which combinations of parameters these distributions can 
probe. We consider it instructive to understand which of these many possibilities have 
been realized in benchmark models, and hope that the present work can help guide future 
studies. 

6 Appendix 

r decays to a single charged particle (plus any number of neutral particles) constitute 84.7% 
of all T decays. Included in this 84.7% is the leptonic branching fraction, 35.2% of all r 
decays, into the final states ez/gZ/^ and ^v^Vt- The remaining 49.5% of all r decays are 
hadronic one-prong rs. Three-prong r's contribute another 15.2%. All branching fractions 
are taken from the PDG ^6] . 

Hadronic one-prong rs offer the most sensitive polarimeter for polarization analyses 
using invariant mass distributions and are the focus of the present work. The dominant 
contributions to the hadronic one-prong decay mode come from 

• — > VTX^: the branching fraction into this mode is 11.1%. 

• r~ — * vp~ VTi~Ti^: the branching fraction into this mode is 25.4%. 

• ^ va^ — s> UTi^Ti^Ti^: the branching fraction into this mode is more difficult 
to estimate. By isospin, this mode is related to the three-prong decay mode — >■ 
uai —>■ VK^-K^-K^ . The total branching fraction into this final state, with no additional 
neutrals, is approximately 9%. We take this to be our estimate of the total branching 
fraction of va^ uti^ti^tt^. 

These three decay modes contribute a branching fraction of 45% to all tau decays, leaving a 
total branching fraction of 5% into one-prong taus from other modes, e.g. modes including 
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kaons and continuum modes with no identifiable intermediate hadronic resonance. We 
model the hadronic one-prong rs as the sum of the contributions from the three decay 
modes listed above. This incorporates 90% of the contributions to the one-prong mode, 
and provides a good approximation to the full hadronic one-prong rs. Incorporating the 
intermediate vector mesons p, Oi is important as the thresholds associated with the masses 
and widths of these particles lead to notable features in the visible hadronic energy spectra 
of polarized rs. 

A full experimental study would be done with TAUOLA [29], which incorporates all 
measured contributions to the r decay modes. 

The matrix elements for the process t —>■ utt depend on the polarization ± of the parent 
r. They are given by (we drop overall constants) 

g 

M+ = COS- (46) 
2 

9 

M- = sin-, (47) 

where 6 is the angle between the pion and the r axis of polarization in the r rest frame. 
Squaring the matrix element and integrating over final state phase space, we obtain 

^ dV 1 ^ 

r a COS 2 

where "Pr = ± is the polarization of the r. These angular distributions can be simply 
understood by appealing to angular momentum conservation: the left-handed neutrino 
must carry off the total angular momentum of the r, and is therefore emitted forward for 
initial left-handed rs, and backward for initial right-handed rs. 

The rs produced at colhders are typically highly boosted, and the variable of interest is 
not 6 but rather z, the fraction of the r (lab frame) energy which is carried by the hadronic 
daughters. 

In the collinear approximation this quantity is invariant under boosts along the r's direction 
of motion. In terms of the (unknown) boost (3 between the r rest frame and the lab frame, 

Pil- {ml/ml))- ^^'^ 

We will henceforth work in the collinear limit, (3^1. Taking this limit and dropping the 
factor m^/rn^, we obtain for the pion spectra 

P±(^) = l±(22;-1). 
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The decay modes r 



up 



T 



uai are more complicated, due to the different 



contributions from the longitudinal and transverse polarizations of the vector mesons, and 
the finite widths of the intermediate states. 

The matrix elements governing r decay to uv^ depend on the polarization ± of the 
parent r as well as the polarization T, L of the vector meson v^. These matrix elements 
are (we again drop overall constants independent of m^,) 



where 9 is the angle between the vector meson and the r axis of polarization in the r 
rest frame. Note the longitudinal polarization L, which carries zero angular momentum 
along the tau polarization axis in the tau rest frame, contributes like the scalar pion, while 
the transverse polarization T contributes oppositely. (Note also that angular momentum 
conservation allows only one of the two transverse polarizations to contribute to the decay.) 
In the narrow width approximation, the net distribution of events in cos 6' is then of the 
form [To] 



In other words, the two different vector meson polarizations make contributions to the 
r helicity-dependent portion of the amplitude which tend to cancel. The factor (m^ — 
2m'l)/{m^ + 2m1) is approximately 0.45 for the p, and only 0.03 for the ai. Thus, with no 
discrimination between the separate contributions of longitudinal and transverse ais, the 
OiS do not yield much sensitivity to r polarization. 

In order to obtain a realistic spectrum of hadronic daughter energies, the finite widths 
of the intermediate vector mesons p, ai must be taken into account. Following [12], we use 
a Breit-Wigner distribution for the vector mesons, incorporating a running width r^,(m^), 
as we now detail. 

In general the contribution to the decay distribution from the intermediate vector meson 
V takes the form 



Here 'P^'^{ml) is the vector meson propagator. It is convenient to write the total matrix 
element for this process as a sum of contributions from the longitudinal and transverse 
polarizations of the intermediate vector meson v. After integrating over the n-pion phase 
space, the interference between the vector meson polarizations vanishes and the decay rate 




(49) 




oc / cin2(r ^ uv) dml dUniv nn) \Mf,{T uv)V^''{ml)Mu{v nn)\'^. 
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given above greatly simplifies. It is natural to define the lineshape (or vector meson decay 
rate) 

Qviml) = J dlln{v rni) \M{v nn)\'^, 

which does not depend on polarization. By the optical theorem the lineshape is proportional 
to the running width of the vector meson, 

m„r„(m^) oc Qviml). 

The decay distribution can now be written 

dV (x ( dIi2{T pv) dml (\Ml{,t uvl)]'^ + {Mrir uvt)]'^) D^{ml)g^{ml), 



where Dy{ml) is the Breit-Wigner with a running width, 



K - ^o) + {mo^oY X 



For the p, we take mo and Tq to be 



(50) 
(51) 



mo,p = 770 MeV, 



To,p = 150 MeV. 



For the ai, we take mo and Fq to be 

mo,a = 1.22 GeV, Fo,a = 420 MeV. 

For the decay p 2tt, the simple chiral Lagrangian interaction 

'Cint oc p^{nidf,Ti2 - vr29/,7ri) 
suffices for our purposes and gives the line shape 

fm^ - 4m!)^/2 



const. X 



m 



(52) 



The parameterization of ai decay, which proceeds dominantly through ai ^ pn ^ Svr, 
is more involved due to the multiple possible parameterizations of intermediate hadronic 
resonances in the chiral Lagrangian describing ai dynamics. Again following [12], we use a 
parameterization of the running width galm^) due to 



ga[m 



4.1(m^ - 9mlY (1 - 3.3(m^ - 9m^) + 5.8(m^ - 9m^)^) m^ < (m^ + m^) 



m 



(1.623 + 



10.38 _ 9^ 



m^ > (m„ + m^r)^ 



(53) 
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All masses are in units of GeV. This particular parameterization follows from a model of 
ai decay which goes through intermediate finite-width p alone, with no contribution from 
other resonances such as the radial excitations p',p"- The coefficients are obtained from a 
fit to experimental data [30] . 

Plugging in the matrix elements (H9l) for the initial decay t ^ vv and doing some 
simplification, we have 



dV (x I d cos 6 dm,, (m^ — m. 



2\2 



mi \mi , 



Dv{ml)g^{ml). 



To change variables from cos to 2, use (HSjl . with m^ in place of m^. With this substitution, 
we have 



dV oc 



1 nzm^ 

dz / dm. 

J (nm-^)'^ 



m, 



m 



2 + 



mz 



mt 



Vr {2Z - I] 



ml 



mt 



mz 



mt 



X D^{ml)g^{ml). 

Finally then, we can write the contribution from the vector mesons as 



1 dV 

r dz 



const. X / dm^ 



1 - 



mz 



2 + 



mz 



mz 



mt 



mt 



+ (2z-l)--f 



mt 



mz 



mt 



-2 



X D^{ml) g^{ml). 



(54) 



Using the lineshapes (!52|) , (!53l) and carrying out the integral over numerically gives the 
hadronic energy spectra PA(z), P^)iz) which are plotted in Figure [T71 
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